Hydroxyl vacancies in single-walled aluminosilicate and aluminogermanate nanotubes 
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We report the first theoretical study of hydroxyl vacancies in aluminosilicate and aluminoger- 
manate single-walled metal-oxide nanotubes. The defects are modeled on both sides of the tube 
walls and lead to occupied and empty states in the band gap which are highly localized both in 
energy and in real space. We find different magnetization states depending on both the chemical 
composition and the specific side with respect to the tube cavity. The defect-induced perturbations 
to the pristine electronic structure are related to the electrostatic polarization across the tube walls 
and the ensuing change in Br0nsted acid-base reactivity. Finally, the capacity to counterbalance 
local charge accumulations, a characteristic feature of these systems, is discussed in view of their 
potential application as insulating coatings for one-dimensional conducting nanodevices. 

PACS numbers: 31.10.+z,31.15.A-,31.15.E-,31.70.Ks 



Since the discovery of carbon nanotubes by Ijima in 
1991 [l[, the field of carbon-based and inorganic nan- 
otube materials has grown considerably due to their 
potential applications in electronics, photonics, chemi- 
cal separation, (photo)catalysis and biotechnology Q. 
While carbon-nanotubes and their inorganic counter- 
parts such as BN, WS2 and M0S2 are routinely produced 
by electric-arc discharge, chemical vapor deposition, or 
laser ablation new low-temperature solution-phase 
chemistry routes have recently been proposed for semi- 
conducting single-walled aluminosilicate (Al-Si) and alu- 
minogermanate (Al-Ge) nanotubes (see Fig. 1 and 0]). 
The reported extremely high experimental control in 
terms of dimensions and monodispersity of these nanos- 
tructures together with the potentially huge range 
of tunable properties via chemical functionalization and 
substitutional doping make both Al-Ge and Al-Si at- 
tractive candidates as large-storage chemical devices Q , 
artificial ion-channel systems Q , and insulating coatings 
for conducting cores, which could make it possible to 
enforce one-dimensional anisotropic conductivity at the 
nanoscale Q. Al-Si and Al-Ge nanotubes are struc- 
turally analogous to the naturally occurring hydrous- 
aluminosilicate Imogolite . Its walls consist of a single 
layer of octahedrally coordinated aluminium hidroxide 
(gibbsite), with tetrahedral silanol (Si-OH) groups at- 
tached at the inner side of the tube (see Fig. 1). From 
a compositional point of view, the only difference be- 
tween Al-Si and Al-Ge tubes is the substitution of silanol 
groups with germanol (Ge-OH) moieties. The resulting 
chemical formula of the unit cell is (A^SiC^H^N and 
(AUiGeOyH^N for the Al-Si and Al-Ge tubes, respec- 
tively. N refers to the number of radially inequivalent 
aluminum atoms along the nanotube circumference, nec- 
essarily an even number for symmetry reasons 0- Re- 



gardless of the specific route adopted for their synthe- 
sis, both Al-Si and Al-Ge tubes are achiral and analo- 
gous to zig-zag (n, 0) semiconducting carbon-based nan- 
otubes 0. Following a recent Density Functional The- 
ory (DFT) study of Imogolite-based nanotubes Q , which 
opened up the possibility of first-principle studies for 
these systems, as well as previous simulations of neutral 
paramagnetic defects in zeolites [Io[ we address the ef- 
fects of the simplest defects, i.e. neutral hydroxyl (-OH) 
vacancies, on Al-Ge and Al-Si in terms of both electronic 
structure changes and ensuing modifications in the global 
insulating properties. 

DFT simulations were performed within a plane-wave 
ultrasoft pseudopotential framework as currently imple- 
mented in the VASP code [ll|]. On the basis ofprevious 
successful results for Al-Si(Ge) based systems [9j, alumi- 
nosilicate mesostructures |l2[ and aluminium oxide thin 
films [H, exchange and correlation were approximated 
at the semi-local generalized gradient-corrected PW91- 
GGA level Q. Hybrid B3LYP [H periodic calcula- 



tions were carried out with the Gaussian p rogr am [lq |. 
adopting a CEP-121G triple-zeta basis set [17]. In the 
spirit of (and at odds with we considered and 
optimized (see Supporting Information, SI) the repeat 
unit cell (c in Fig. 1) for Al-Si (Al-Ge) systems made 
up of N — 24 (N = 36) gibbsite units. The optimized 
c value of 8.68±0.01 A (8.78±0.01 A) is slightly larger 
than experimental X-ray diffraction (XRD) data for Al- 
Si (Al-Ge) thin films i.e. 8.51 A [|. The deviations 
are due to the well-known overestimation of aluminosil- 
icate bond len gths by the semi-local density functional 
approximation 12, 18, T^. We also note that the agree- 
ment between the calculated values of 8.68 A for N = 24 
considered here and the reported 8.62 A for N=20 (equiv- 
alent to N u =10 in '{|) confirms that the size of the unit 




FIG. 1: (Color online) Optimized geometrical structure of 
(Al 2 Si0 7 H 4 )24 (left) and (Al 2 Ge0 7 H 4 )36 (right) based nan- 
otubes. The single-wall structural motif (bottom) is displayed 
together with the zig-zag periodic unit of size c, cLai-ai along 
the nanotube axis and circumference. O: red, H: gray, Al: 
green, Si: yellow, Ge: cyan. 
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cell depends only weakly on N [9j . The optimized diam- 
eter of 23.15±0.05 (33.1±0.04) A is in good agreement 
with the reported experimental value of ^22 (~33) A 
for Al-Si (Al-Ge) By independently optimizing the 
bi-dimensional (2D) analogues of Al-Si and Al-Ge (see 
SI), the bending energy was evaluated to be 0.45 eV/N 
and 0.19 eV/N for Al-Si and Al-Ge tubes respectively. 
These results further confirm that the tube bending is a 
consequence of minimizing the strain from the mismatch 
between the stronger Si(Ge)-0 and weaker Al-0 bonds 
(see SI and 0,01) in that the more bent Al-Si is more sta- 
bilized than Al-Ge by the bending process. We also note 
that on the basis of the calculated smaller Al-Ge forma- 
tion energy, it is tempting to link these results with the 
experimental evidence of longer (~100 nm) Al-Si tubes 
and shorter (~20 nm) Al-Ge tubes 0]. 

The calculated Density of States (DOS) for Al-Ge, Al- 
Si and the corresponding 2D sheets are reported in Fig. 
2. The calculated band gaps are 4.1 eV and 3.9 eV for Al- 
Si and Al-Ge respectively. Besides being in close agree- 
ment with the experimentally reported value of 3.6 eV 
for Al-Ge : 3j, these values are in accordance also with 
simulations of aluminosilicate clusters (~4 eV [l2| ) and 
aluminum oxide thin films (~4 eV [HI) obtained using 



FIG. 2: (Color online) Total Density of States (DOS) for 
defect-free Al-Si and Al-Ge, their 2D analogues, and in the 
presence of one OH vac both on the outer (OUT) and the in- 
ner (IN) surface of the tubes. Calculated Fermi energies are 
displayed as dotted line with the same color labeling as for 
the DOS. 2D and band gap defect states (filled) have been 
increased by a factor 10 for clarity. 

the same PW91 functional. Finally, we note that all data 
match also the calculated band gaps recently reported for 
(N=20) Al-Si, Al-Ge systems [9J. When considering the 
2D-sheets, the physical bending of the nanotubes is found 
to affect only Al-Si but not Al-Ge. In fact, while we do 
not find any significant change in the calculated band 
gap between nanotube (3.9 eV) and 2D-sheet (4.0 eV) 
for the less bent Al-Ge tube, for the Al-Si tube, due to 
the smaller radius of curvature, the calculated band gap 
is 0.4 eV smaller than for the corresponding 2D-sheet. 
For completeness we considered the Al-Si (Al-Ge) 2D- 
sheets also by adopting the hybrid functional B3LYP [l5[ 
which is reported to account correctly for the chemical 
reactivity of aluminosilicate-based materials In line 
with previous HF-based calculations for zeolites, which 
produced band gaps considerably larger than PW91 re- 
sults the B3LYP band gaps are wider than both 
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the PW91 values and the experimental optical gap of 3.6 
eV Specifically, within the adopted triple-zeta basis 
(see SI), we found 9.3 eV and 7.42 eV for Al-Si and Al-Ge 
2D-sheets respectively. On these grounds it is reasonable 
to expect that B3LYP will also overestimate the band 
gap for the complete nanotubes. 
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FIG. 3: Band-decomposed charge density for defect-free Al- 
Si (a) and Al-Ge (d), and in the presence of one OH vac 
both inside (b,e) and outside (c,f) the nanotube cavity. The 
VB (green, l(r 6 eA- 3 ) and CB (red, 5 x 10 _7 eA- 3 ) densi- 
ties have been integrated over 0.5 eV from the band onset. 
Occupied and empty band-gap defect states are displayed 
(5 x l(T 7 eA- 3 ) in cyan and pink, respectively. 



Fig. 3 shows the band-decomposed charge densities for 
both the valence band (VB) and conduction band (CB) 
edges, calculated by considering electronic states within 
0.5 eV from the band onsets. Interestingly, it emerges 
that both Al-Si and Al-Ge are characterized by a neat 
separation in real space of VB and CB. The VB edge 
is in fact localized inside the nanotube cavity, while the 
CB edge faces the outer side of the nanotube. In line 
with [9J] , these findings confirm a strong radial anisotropic 
electron affinity with an ensuing enhanced Br0nsted acid- 
ity (basicity) for the outer (inner) tube surfaces as sug- 
gested in [22| (see SI for further details). In order to 



assess the effect of OH vacancies on the global electronic 
structure of the tube, one hydroxyl fragment was elimi- 
nated both on the inner, Al — Si(Ge) — OHj nvac , and on 
the outer, Al — Si(Ge) — OH outV ac, side. In anal ogy with 



O vacancies on other metal-oxide substrates [23j, the 
presence of one OH vac is found to introduce electronic 
states in the pristine band gap (see Fig. 2). One OHi nvac 
is found to create both occupied and unoccupied defect 
states, which reduce the actual band-gap to ~1.8 eV and 
eV for both Al-Si and Al-Ge. This effect thus en- 
hances the conductivity for the global systems substan- 
tially. Both defects are magnetic, generating a doublet 
spin-state. In line with previous results for paramagnetic 
defects in zeolites [lfj, both for Al-Ge and Al-Si, the oc- 
cupied (unoccupied) defect states (see Fig. 3 and SI) 
are highly localized around the undercoordinated Si (Go) 
atom. Interestingly, while one OH vac on the outer wall 
of Al-Si is also found to create a localized paramagnetic 
(doublet) spin-state and to reduce the actual band gap 
(0.8 eV), the same defect in the Al-Ge analogue forms 
a state just at the CB onset which pins the Fermi level 
there. Consequently, despite the odd total electron count 
in the system, our DFT energy is minimized by a non- 
magnetic solution with equal occupancy of both up and 
down spins [2l[ (see Fig. 3 and SI). We cannot be sure 
from our calculations whether this behaviour would also 
be found for a truly isolated defect, or is an artefact of 
the periodic defect structure along the tube imposed by 
our boundary conditions, but we would expect the bind- 
ing energy of any localized state to be extremely small. 
The calculated energies for both OH outV ac occupied states 
(Fig. 2) are very close, suggesting that the reduced (yet 
still overestimated in PW91 Q) band gap of Al-Ge is the 
cause of this metallization. 

It has been suggested that surface charge properties 
and chemical reactivity of Imogolite depend on a deli- 
cate (pH-dcpendent) balance between local geometrical 
deformations associated with the tube bending and posi- 
tive (ne gati ve) charge accumulation on the outer (inner) 
surface [22|. Charge separation across the tube walls is 
expected to play a fundamental role also for the suggested 
use of Al-Si (Al-Ge) as an insulating coating around con- 
ducting cores 0, Q. In this respect, it is of utmost 
importance to assess the changes induced by OH vac in 
terms of both local and global polarizations across the 
tube section. Following 24(, the change (AV) in the 



(microscopically) averaged electrostatic potential across 
the polarized interface can be related to the electrostatic 
dipole across the interface itself (p) as AV = e/i/eo. On 
this basis, it is straightforward to evaluate the actual 
dipole moment across the tube walls from AV values 
calculated in DFT. From the shift in the electrostatic 
potential (averaged both angularly and along the tube 
axis, see Fig. 4), Al-Ge is calculated to possess a larger 
dipole across the wall (0.676 Debye or roughly one third 
of H 2 0, 1.85 D) with respect to Al-Si (0.591 D). In line 
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FIG. 4: (Color online) Top: radial distribution of the averaged 
electrostatic potential V (see text for details of averaging). 
Bottom: radial V shifts as a function of the considered OH vac 
defects. Vertical arrows highlight the electrostatic potential 
shifts AV across the tube walls. 

with [23] (and at odds with 0) the calculated dipole di- 
rections agree with the positive (negative) charge accu- 
mulation at the outer (inner) surface. Despite its larger 
dipole, Al-Ge is characterized by lower inner and outer 
plateaus of V with respect to Al-Si. This suggests that 
the net polarization comes from a delicate balance in- 
volving the total electronic distribution across the whole 
interface, rather than the atomic charges alone. Owing 
to the angular symmetry of the tubes, both Al-Ge and 
Al-Si do not possess any permanent global dipole mo- 
ment. Averaging [25| the electrostatic potential over the 
circular sector affected by the defect states (see Fig. 3), 
it turns out that the introduction of one OH outV ac with 
the ensuing electronic charge accumulation on outer wall 
slightly lowers the corresponding dipole by 0.021 D (0.028 
D) for Al-Si(Ge). The results are less intuitive (at least 
with respect to the defect induced charge accumulation 
in Fig. 3) for OH invac defects, where the correspond- 
ing dipole is increased (decreased) by 0.014 D (0.008 D) 
for Al-Si(Ge), further suggesting that the total defect- 
induced redistribution of electron charge must be taken 
into account even for the highly localized perturbations 
found here. Interestingly, the lack of global polarization 



is maintained for both Al-Ge and Al-Si even in the pres- 
ence of one OH vac . These findings suggest a capacity of 
both Al-Ge and Al-Si to counterbalance intrinsically lo- 
cal charge accumulations and to damp any ensuing global 
polarization, which could be highly relevant for potential 
applications as insulating components in one-dimensional 
conducting-semi-conducting hybrid nanodevices. 
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